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We investigate the stability of a quadratic band-crossing point (QBCP) in 2D fermionic systems. 
At the non-interacting level, we show that a QBCP exists and is topologically stable for a Berry 
flux ±27r, if the point symmetry group has either fourfold or sixfold rotational symmetries. This 
putative topologically stable free-fermion QBCP is marginally unstable to arbitrarily weak short- 
range repulsive interactions. We consider both spinless and spin-1/2 fermions. Four possible ordered 
states result: a quantum anomalous Hall phase, a quantum spin Hall phase, a nematic phase, and 
a nematic-spin-nematic phase. 
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Introduction — In multi-band fermionic systems, a 
band-crossing point (BCP) is a point in the Brillouin 
zone where two bands cross. As the chemical poten- 
tial reaches a BCP, the Fermi surface shrinks to a point 
and new phenomena, not described by a Fermi liquid, re- 
sult. The simplest and best studied is the case of a linear 
band crossing, whose low energy physics is described by 
a Dirac fermion. Dirac fermions are a good description of 
the low energy states of nodal superconductors, graphene 
and zero gap semiconductors. 

In general, a Dirac point in a band structure is ro- 
bust with respect to small changes in the effective poten- 
tial which preserve the symmetries of the crystal, as has 
been extensively shown in various contexts [l|, [3, H, 
Moreover, short-range electron-electron interactions are 
perturbatively irrelevant for space dimension d > 1. 
Thus, there is a stable phase with free gapless Dirac 
fermions which becomes unstable above a critical interac- 
tion strength corresponding to a quantum critical point 
beyond which lie phases with spontaneously broken space 
or point group symmetries and/or broken time-reversal 
invariance 0, @, [3|- 

In this letter we consider a system with a quadratic 
band crossing point (QBCP) somewhere in its 2D Bril- 
louin zone. This problem has not been discussed in 
depth, and only a few aspects have been analyzed. The 
perturbative stability of a QBCP was studied for 2D 
noninteracting systems with symmetry in Ref. Q. 
For interacting fermions, it was noted in Ref. 0] that a 
QBCP in 2D has instabilities, for arbitrarily weak inter- 
actions, leading to the spontaneous breaking of rotational 
symmetry (nematic phase) or time-reversal invariance, 
but its consequences were not explored in depth. 

We begin by analyzing the general symmetry principles 
that protect a QBCP in lattice models of noninteracting 
fermions. We find that QBCPs are protected by time- 
reversal symmetry and C4 or Cg rotational symmetry. 
Explicit examples of lattice models with both symme- 



tries are presented. We show that short-range repulsive 
interactions are marginally relevant in the renormaliza- 
tion group (RG) sense. The symmetry breaking phases 
and phase transitions, both quantum and thermal, are 
investigated by an RG analysis and mean-field approxi- 
mations, both presumably reliable at weak-coupling. 

We determined the structure of the phase diagrams for 
both spinless and spin-1/2 fermions. In the spinless case 
the leading weak-coupling instability is to a gapped phase 
with broken time-reversal invariance, a quantum anoma- 
lous Hall (QAH) effect, and topologically protected edge 
states. For stronger interactions, there is a subsequent 
transition to a nematic Dirac phase and an intermedi- 
ate phase with QAH-nematic coexistence. For spin-1/2 
fermions the phase diagram is more complex: in addition 
to spin singlet QAH and nematic phases, there are also 
a spin triplet quantum spin Hall (QSH) phase[l3] and a 
nematic-spin-nematic (NSN) phase [HI, El, [13] ■ 

Quadratic band-crossing point — A BCP carries quan- 
tized Berry flux [l^ as required by time-reversal symme- 
try: — i dk- (-(/^(k)! Vk|^/'(k)) = nTT, where r is a contour 
in the momentum space enclosing the BCP, V'(k) is the 
Bloch wave function in a band involved in the band cross- 
ing, and n is an integer. For a Dirac point, the Berry flux 
is ±7r. Instead, the Berry flux at a QBCP is either or 
±27r [§| . The zero flux QBCP is an accidental band cross- 
ing, which can be removed by infinitesimal band mixing 
without breaking any symmetries, but a QBCP with ±27r 
flux is robust and more interesting. 

A natural question to ask is if a QBCP is protected 
by the symmetries of the non-interacting system. In gen- 
eral, there are two ways to remove a QBCP. One way 
is to split it into several Dirac points while preserving 
the total Berry flux. A QBCP with flux 27r, for example, 
can be split into two separate Dirac points each with flux 
TT, or three Dirac points with flux tt and one additional 
Dirac point with —tt. The former case in general breaks 
the point group symmetry leaving, at most, a two- fold 
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rotational symmetry unbroken. The latter case can take 
place while preserving a three-fold rotational symmetry, 
such as the case of bilayer graphene 3, lB| • For a QBCP 
with a fourfold or sixfold symmetry axis, the split into 
Dirac points cannot occur without breaking that sym- 
metry. The alternative is to open a gap by breaking 
time-reversal symmetry or a symmetry that is formally 
similar, such as the combined space and spin symmetry 
whose breaking leads to a spontaneous quantum spin Hall 
state, as discussed below. Thus, for a QBCP (with Berry 
flux ±27r) to be stable without fine tuning, two conditions 
are required: a) the system must be time-reversal invari- 
ant and b) the QBCP must have C4 or Cg symmetry. 

An example of a QBCP in 2D with C4 symmetry can 
be found in the checkerboard lattice Fig. [T(a)| This 
lattice can be regarded as the 2D projection of a 3D py- 
rochlore lattice. It is also the oxygen lattice in a Cu02 
plane of the cuprates. With one orbital per site, there 
are two bands crossing at a QBCP at (tt, tt) with a four- 
fold rotational symmetry. At half filling the QBCP is at 
the Fermi level. An example of a QBCP with Cq symme- 
try is a tight-binding model on a Kagome lattice, which 
has three bands. The middle band touches the bottom 
band at (0,0), resulting in a QBCP with sixfold rota- 
tional symmetry. It lies at the Fermi level at 1/3 filling. 



charged Fermi fields, tpi and ■02, whose Hamiltonian is 
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FIG. 1: (Color online) (a) A checkerboard lattice and (b) a 
Kagome lattice. The arrows represent currents in a sponta- 
neously generated QAH state that breaks the time-reversal 
symmetry. See text for details. 

In the presence of weak interactions, a BCP may be- 
come unstable if interactions are relevant in the RG sense. 
In 2D, a QBCP has a finite one-particle DOS, which im- 
plies that short-range interactions are marginal at tree 
level. We will show below, that at a 2D QBCP a short- 
range repulsive interaction is marginally relevant, and 
destabilizes this free fermion fixed point in weak coupling, 
leading to a state which spontaneously breaks one of the 
symmetries that otherwise would protect the QBCP. 

General model — We first formulate the theory of pos- 
sible symmetry breaking phases in a general way. We 
begin with the spinless-fermion case. Near a QBCP, in 
the low energy regime we have two species of interacting 



H= dv 



*t(r)Ho*(r) + y^{(r)V;i(r),A|(r)^A2(r) , (1) 



where = {tp\,^\), '4' is its conjugate, and V is the 
coupling constant of the interaction. 

The band structure near the QBCP is obtained by di- 
agonalizing 2x2 Hermitian matrix 7io(k) for all Bloch- 
wave vectors in the neighborhood of the band crossing 
point, |k| <C 1. Quite generally we can choose the iden- 
tity matrix / and the two real Pauli matrices cr^ and 
as a basis [Tt'I and write Ho(k) as 18 1 



7i:o(k) = dil + d.^a'' + d^a^ 



(2) 



where di = ti{k1 + ky), d^ 
tz{k'^ — ky). The d-wave symmetry of dx and dz dis- 
tinguishes a QBCP from a Dirac point in which their 
counterparts have a p-wave symmetry. It is this d-wave 
nature that gives rise to the ±27r Berry phase of a QBCP. 
For a QBCP with a Cq rotational symmetry, \tx\ = \tz\- 
If the system has particle- hole symmetry, i/ = 0. The 
condition \ti\ < \tx\ and < \tz\ is required to ensure 
that away from the QBCP, one of the bands lies above 
the degenerate point and the other band lies below. 

At y = 0, in the model of Eq. ([T]), the fermions have 
a finite DOS but do not have a Fermi surface. They 
have a dynamic critical exponent z = 2, and an effec- 
tive dimension dcff = d + z = 4 [foj . 4* has dimension 
one ['^] — 1, and the only local four- fermion operator 
allowed is marginal since 4[^] = dcs- There is a single 
dimensionless coupling constant g = V/\tx\. 

This system is similar to = 1 spinless fermions, a 
system with two Fermi points and dynamic critical expo- 
nent z = 1. In the ID case the Fermi field has scaling di- 
mension 1/2, so there is only one interaction, four- Fermi 
backscattering, which is potentially important. However, 
due to a cancellation between the Cooper channel and the 
bubble term in ID, the interaction is exactly marginal to 
all orders in perturbation theory [2Q||, which is the origin 
of Luttinger liquid behavior in ID. In contrast, no similar 
cancellation occurs for fermions in 2D with z ~ 2. Al- 
though the 4- Fermi interaction is superficially marginal it 
is actually marginally relevant. We find, that to one-loop 
order, the RG beta function for g ~ V/\tx\ is 
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where 



is a momentum 



[vi-{tz/txr), I is 

rescaling k — > fce^', and K{x) is the complete elliptic 
integral. For \tx \ = \tz\, i-e. a QBCP with Cg symmetry, 
a = (47r)~^ Hence, Eq.Q imphes that for g > 

the effective coupling constant flows to strong coupling. 

To explore the consequences of this instability, we in- 
vestigated, in a mean-field level, possible orderings of bi- 
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linear order parameters: 

(*^(r)ay*(r)), 



(4) 



$ is the order parameter of a time-reversal symmetry 
breaking gapped QAH phase 0, 01 • This phase has a 
zero-field quantized Hall conductivity axy — e'^/h. Qi 
and Q2 describe the nematic phases in which the C4 or 
Cg rotational symmetry is broken down to C2 by splitting 
the QBCP into two Dirac points located along the direc- 
tion of one of the main axes (Qi), or along a diagonal 
iQ2)- The nematic phase is an anisotropic semimetal. 
Unlike in graphene, where the two Dirac points have 
Berry fluxes tt and — vr, in the nematic phase both Dirac 
points have the same Berry flux. There is also a phase 
in which nematic {Qi ^ or Q2 ^ 0) and QAH or- 
ders ($ ^ 0) coexist, an insulating analog of the metallic 
time-reversal breaking nematic /3 phases of Ref. ^ . 

Since there is only one coupling constant {V) in Eq.dl]), 
the weak-coupling ordering tendencies are determined by 
the logarithmically divergent normal state susceptibilities 
X$ (QAH order) and XQi and XQ2 (nematic order). For 
general and t^, they satisfy x* = XQi + XQ2- Hence, 
X* > XQi (* = Ij 2), so the leading weak coupling insta- 
bility is to the (gapped) QAH state. 

The mean-field Hamiltonian is 



MF 



V 



*(r) 

(5) 



By minimizing the ground state energy of -ffiviF we find 
that at weak coupling the ground state is indeed the QAH 
phase, with a gap A ~ Aex-p{—2/ag) (A is a cutoff) and a 
mean-field critical temperature Tc ~ A, consistent with 
the scaling predicted by the RG. A 3D example of the 
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QAH at finite coupling is discussed in Ref. 

Mean-field theory also predicts nematic phases 
provided that irrelevant operators, such as 
Jdvdv'j:^^,,Uir - r')^l(r)V.(r)Vj(r')V'.(r'), are 
also included. The nematic phase Q\ is energetically 
favored at small F > and t/ < if |C//^| is large 
enough. As |J7/V^| is reduced, the nematic phase gives 
way to the QAH phase (and to a mixed phase). 

Lattice models — We consider the following minimum 
model on a checkerboard lattice with a QBCP. 



H 



E 



(6) 



where tij is the hopping amplitude between sites i and 
j and > is the nearest-neighbor repulsion. Here, 
tij — t, t', t" , respectively for nearest neighbors, and next- 
nearest neighbors connected (or not) by a diagonal bond. 
Fig. 1(a) There are two sublattices A (red) and B 



(blue). The fermion spinor is = (c'^jC'b)- param- 
eters of the free fermion Hamiltonian [Eq. ([2])] are dj = 
— [f + t") [cos kx + cos ky), dx — — 4icos^cos^, and 

4 = -{f -t"){cosky-coskx). The QBCP is M= (7r,7r), 
at the corner of the Brillouin zone. The parameters of the 
continuum Hamiltonian (near the QBCP) of Eq. ^ are 
ti = {t' + t")/2, tx = t/2, t^ = (f -t")/2. The order pa- 
rameters are Qi = jJ2si^A,i^A,i - c'BA+s^B,i+s) ("site 
nematic"), Q2 = ^J2s ^s'R.e{c\ -CB,i+s) ("bond ne- 
matic"), and $ = lJ2s ^sl'm-{c\ ^CB,i+s) (QAH), where 

5 = ±x/2 ± y/2 connects nearest neighbors. Ds — ±1, 
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FIG. 2: (Color online) Mean-field T — V phase diagram of a 
half-filled checkerboard lattice with t' /t = 0.5, t"/t = -0.2. 
Nematic and QAH orders coexist in the shaded area. Thick 
(thin) fines are first (second) order transitions. A is a bicritical 
point and i? is a critical-end point. 



A mean- field theory analysis 23[ (similar to Eq.®) of 
the lattice model, Eq.®, yields the T — V phase diagram 
of Fig. H (The details depend on t' /t and t" /t.) A QAH 
phase is found for V small and below a critical temper- 
ature. This phase has a zero-field quantized Hall con- 
ductivity e^/h, and that the quasiparticle spectrum has 
topologically protected chiral edge states, as predicted by 
general considerations 2j| . A site- nematic phase is found 
ior V ^ \t' — t"\, while the bond-nematic is not favored. 
For \t'\ > \t"\ and \t"\/\t'\ 1, there is a direct nematic- 
QAH first order transition. If ~ 1, there is also a 

coexisting QAH-|-nematic phase. For other values, there 
are a direct first-order transition and a coexisting phase. 

Near a QBCP, a next-nearest-neighbor attraction 
V„nn < [25| generates the (irrelevant) interaction U. 
For \Vnnn/V\ 1, a site-nematic phase (Qi) becomes 
stable, even at weak coupling. A site-nematic phase is 
stabilized at large enough V > {Vnnn = 0). In the 
strong coupling limit, at this filling, this system is known 
to be in an insulating site-nematic phase (26l | . 

Spin-1/2 fermions — Let us consider briefly the case of 
a system of spin-1/2 fermions at a QBCP. Details will be 
given in Ref. [1^. We consider the four shortest-range 
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interactions on a lattice: a) an on-site repulsive Hubbard 
U, b) a nearest neighbor repulsion V, c) a nearest neigh- 
bor exchange interaction J, and d) a pair-hopping term 
W. In addition to spin-singlet order parameters {c.f. Eq 
([5])), there are also spin-triplet order parameters: 

Q{ = (*t(r)(r®a.)*(r)), = (*t(r)(r ® a,)*(r)), 
5^= (*t(r)(f 0/)*(r)), $* = (^'t(i.)(f ^cry)*(r)}. (7) 

where r are the three Pauli matrices. Here, 5* is the 
spin density. For Q\ ^ the QBCP splits into four 
Dirac points displaced along the main axes. This state 
has reversed spin polarization along x and y axes. For 
a QBCP with C4 symmetry, the charge sector is still 
C4 invariant, but the spin sector changes sign under a 
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rotation by tt/2. Thus, Q\ 7^ is a NSN state 
Q2 7^ describes NSN order along the diagonals. 

A state with ^ is a QSH phase 0, [3, [H, [lil with 
helical edge states 3^ 31|. In this phase, the two spin 
components have opposite Hall conductivity. The QSH 
topological insulator has a quantized spin Hall conduc- 
tivity, and it is a spin triplet version of the QAH phase. 
It is an insulating analog of the /3-phase of Refs. IJ, 33 |. 

A scaling analysis similar to the spinless case finds six 
marginally relevant operators, associated with the four 
nematic order parameters (singlet and triplet) and the 
QAH and QSH orders. We have obtained the phase di- 
agram in mean-field theory for spin- 1/2 fcrmions on a 
checkerboard lattice at a QBCP with U > V > at 
T = 0[33]. For simplicity we take t, t' and t" such that 
tx ~ ^z, and W = Q. Once again all the susceptibil- 
ities are logarithmically divergent with x$ = X*: and 
Xq. — XQi (* — 1:2), so the interactions are marginally 
relevant. For U > 2V > 0, the system is in the NSN 
phase at low temperatures. For 2V > U > we find the 
QAH phase for J > 0, and the QSH phase for J < 0. 
Both gapped phases are topological insulators. The gaps 
and critical temperatures obey a scaling law similar to 
the spinless case [2^. For U ^ 00 (and J = 0), there is 
a NSN state, while U — > 00 stabilizes a nematic phase. 

Using RG methods and mean-field theory we showed 
that a system of interacting fermions, with or without 
spin, at a QBCP have topological insulating QAH or 
QSH phases, at arbitrarily weak short-range repulsive 
interactions. These perturbatively accessible topologi- 
cal insulating phases are due to spontaneous symmetry 
breaking, described by order parameters, and are not due 
to spin orbit effects in the band structure. At intermedi- 
ate coupling we also find nematic (and coexisting) phases. 
Using large N methods and a 2 -I- e expansion, we infer 
the existence of these phases in 3D (similar to those of 
Ref.ji^l), but at a finite critical coupling [isj. 
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